Recoil Order Chiral Corrections to Baryon Octet Axial Currents 
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Abstract 



We calculate chiral corrections to the octet axial currents through 0(p 3 ) 
^-v \ using baryon chiral perturbation theory (BCPT). The relativistic BCPT 

framework allows one to sum an infinite series of recoil corrections at a given 
order in the chiral expansion. We also include SU(3)-breaking operators oc- 
^^ ■ curing at 0(p 2 ) not previously considered. We determine the corresponding 

low-energy constants (LEC's) from hyperon semileptonic decay data using a 
variety of infrared regularization schemes. We find that the chiral expansion 
of the axial currents does not display the proper convergence behavior, regard- 
less of which scheme is chosen. We explore the implications of our analysis 
SZ • for determinations of the strange quark contribution to the nucleon spin, As. 

>■ 

PACS Indices: 12.39.Fe, 13.30.Ce, 14.20. Jn 



I. INTRODUCTION 

In the pseudoscalar meson sector chiral perturbation theory (CPT) provides a consis- 
tent and systematic framework for calculating physical observables. Generally, they can be 
expanded order by order in powers of p/A x , where A x = 4:iiF n , F w = 92.4MeV, and p is the 
typical small momenta or mass of the Goldstone bosons. When CPT is extended to include 
the baryons, a difficulty arises due to the presence of the large baryon mass. One encounters 
terms like m^/A x , which obscure the convergence of the chiral expansion. To overcome this 
difficulty, heavy baryon chiral perturbation theory (HBCPT) was introduced @@,i]- In this 
non-relativistic framework, the baryon mass appears only in vertices as powers of l/m^. 
One thus obtains a consistent expansion in two small parameters, p/m^ and p/A x . This 
approach has been applied successfully to a wide variety of baryon observables. 

Despite its successes, HBCPT comes with its own shortcomings. For example, the I/tun 
corrections are unnaturally large for some observables. In some cases, one requires a large 
number of higher-order terms in l/m^ in order to obtain proper convergence behavior of 
the chiral expansion. From a conceptual standpoint, it was noted in [l[] that HBCPT fails 
to produce the correct the analytical structure near threshold for the nucleon isovector 
electromagnetic form factors. The underlying reason is the non-relativistic treatment of the 
baryon propagators in HBCPT. 

A relativistic formulation of CPT for baryons was recently proposed in Ref. || and 
applied to nucleon electromagnetic form factors in Ref. J7J. This formulation, which we 
denote BCPT (baryon CPT), circumvents the problematic {ttin/ A x ) n terms by splitting a 
given chiral loop integral into an infrared sensitive term, J, and an infrared insensitive, or 
"regular" piece, R. The former contains all the non-analytic contributions uniquely identified 
with chiral loops; the latter contains the power dependence on baryon mass. Since R is also 
analytic in quark masses, its contribution may be completely absorbed into the appropriate 
terms in a chiral Lagrangian. Since the corresponding low-energy constants (LEC's) are 
determined entirely from fits to experimental data, the (mN/A x ) n behavior never appears 
explicitly. Moreover, by retaining the fully relativistic form of the baryon propagators in /, 
one includes all of the recoil corrections to the non-analytic contributions at a given order in 
the chiral expansion. This procedure, known as "infrared regularization" , contrasts with the 
HBCPT approach, where one must explicitly work out the recoil corrections order-by-order 
in 1/rriN- 

The simplifications introduced by BCPT have been explored in the case of a few observ- 
ables. In Ref. 0, for example, it was pointed out that BCPT improves the convergence of 
the chiral expansion of the nucleon electromagnetic form factors as compared to HBCPT. 
Moreover, since BCPT is relativistic, analytical behaviour of the resulting form factors is 
correct. 

In this work we employ BCPT to calculate the one-loop chiral corrections to the axial 
currents of the octet baryons. The leading order chiral corrections to the axial currents are 
of the form m s In m s and were first calculated in |2j . For a subtraction scale of \i = 1 Ge V 
the corrections calculated in [0 are less than 30%. However, the correction due to the wave 
function renormalization was ignored in as pointed out in ||, where the same problem was 
treated with HBCPT formalism. When wave function renormalization and vertex corrections 
are both included, the leading one-loop correction is large M. For example, the fit values of 



the SU{3) couplings at tree level in Ref. |] are D = 0.80 ± 0.14, F = 0.50 ± 0.12. The one- 
loop chiral corrections shifted the best fit to D = 0.56 ± 0.1, F — 0.33 ± 0.06 ||. Later the 
same authors included the intermediate decuplet baryon states in the chiral loops and found 
significant cancellations with the octet contributions ||. While this cancellation suggests 
the importance of including the decuplet for obtaining proper convergence, inclusion of the 
decuplet is not sufficient in the case of some other observables. In the case of octet baryon 
magnetic moments, for example, one must also include the leading recoil-order (l/m^v) 
corrections [§]. 

In this paper, we use BCPT to explore the effect of recoil corrections on the convergence 
of the chiral expansion of the octet axial currents. We also include 0(p 2 ) chiral symmetry 
breaking terms not included in previous analyses. In order to maintain predictive power, 
we truncate the expansion at 0(p 3 ). The number of LEC's appearing at 0(p 4 ) prevents 
one from carrying out a model independent analysis. We also follow Refs. @|| and set 
m u = rrid = in performing numerical fits, although the formulae presented below included 
results for non- vanishing pion mass. We find that the impact of the 0(p 2 ) symmetry breaking 
(SB) operators is noticeable. More importantly, the 0(p 3 ) contributions - corresponding 
entirely to loop-generated recoil corrections - are generally larger than the 0(p 2 ) terms. 
Thus, the chiral expansion of the axial currents does not converge in the manner expected 
when decuplet contributions are integrated out. While the significance of the 0(p 2 ) loop 
corrections was first noted in Ref. 0, our study of the expansion through 0(p 3 ) makes the 
non- convergence of the series abundantly clear. Contrary to one's naive hope, inclusion of 
octet-only recoil order contributions only worsens the convergence properties of the axial 
currents. Whether explicit inclusion of recoil order decuplet contributions remedies this 
situation remains to be seen. 

In a related issue, the definition of the infrared loop contributions contains a degree of 
ambiguity. While the non-analytic quark mass dependence of / is unique, this integral may 
also contain terms analytic in m q . Whether or not one retains these analytic contributions 
explicitly is a matter of convention. The standard practice in HBCPT is to keep only the non- 
analytic loop effects. On the other hand, the authors of Refs. ||,|7j also retain analytic pieces 
of /. We analyze the axial currents using both schemes. In this case, the difference amounts 
only to the treatment of 0(p 2 ) analytic contributions to wavefunction renormalization. The 
corresponding impact on the convergence properties of the chiral expansion is small. We 
cannot, however, determine whether this scheme insensitivity persists to higher order. While 
the integrals / for the axial currents contain a variety of 0(p 4 ) and higher contributions, we 
truncate at 0{p 3 ) for reasons noted earlier. 

II. INFRARED REGULARIZATION 

The motivation and formalism for BCPT are explained extensively in Refs. pM. In- 
terested readers may consult these two references for details. The key feature of BCPT 
is the so-called infrared regularization procedure. Following Refs. |5|,[/] we illustrate using 
the one-loop baryon self energy. The ultraviolet (UV) divergence of the one-loop integral is 
regulated using dimensional regularization. The regulated integral H is then separated into 
the / and R pieces using Feynman parameters: 
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{2ir) d AB Jo J (2n) d [(1 - z)A + zB} 2 

~ H C ~ T) dz I Wy [{i- z) l A + zBY = I + R ' (1) 

with A = M 2 — k 2 — ie, B = m 2 N — (p — k) 2 — ie and M, m^ is the pseudoscalar and nucleon 
mass. The region of parameter integration for the integral I contains z — 0. At the origin, 
the denominator is proportional to A 2 , and thus, is highly infrared sensitive (singular in the 
case of the self energy). As shown in Ref. ||, all of the non-analytic m,j-dependence uniquely 
associated with the loop is contained in I. For the regular part, R, the Feynman parameter 
runs from one to infinity, and the result is analytic. Consequently, its contribution can be 
entirely absorbed into the appropriate operators appearing in the effective Lagrangian. In 
addition, if we expand / in terms of 1/mjv, we recover the HBCPT result at each order. Thus, 
the infrared sensitive part I of the corresponding relativistic diagram is just the sum of the 
leading HBCPT diagram and diagrams with 1/mjv insertions to all orders. In other words, 
BCPT effectively sums the 1/mjv series in HBCPT. 

It is important to note that the inclusion of the full tower of recoil corrections renders 
the chiral counting somewhat ambiguous. Contributions involving recoil effects have the 
generic form 



m x ,,k 



A l 



n k M , (2) 



where mx is a pseudoscalar mass, /x = m x /m N oc ^/FrT^/mx, and /(/i) is a recoil factor. 
For the axial currents, /i fc /(/x) is non-analytic in m q and, therefore, can never be generated 
by terms in the effective Lagrangian. Nevertheless, one may perform a Taylor expansion 
of /(//) in powers of /i 2 about /i 2 = 0. Consequently, the quantity in Eq. (Q) contains an 
infinite series of contributions of successively higher orders in p. The first term in the series 
- obtained by replacing /(//) — > /(0) - is purely of 0(p k+2 ). In the language of HBCPT, 
this first term in the series constitutes its leading contribution in the 1/mjv expansion. In 
what follows, we identify the chiral order of the term in (H) by the order of its leading term 
in the 1/mjy expansion. 

We emphasize that the chiral order of the recoil term in (0) is unambiguous only in the 
heavy baryon limit. Retention of the higher order terms associated with recoil factors is the 
price one must pay for maintaining the analyticity properties of loops implied by relativity 
and crossing symmetry. These properties are lost in HBCPT. It does not appear possible to 
respect the full analytic structure of chiral loops and maintain the standard chiral counting 
procedure simultaneously. Fortunately, in the case of the axial currents, we find that the 
numerical impact of setting /(//) — > /(0) is negligible. In short, it is sufficient to work to 
0(1/ m^) in the heavy baryon expansion in order to ascertain the effects of recoil. 

It was also argued in Ref. that the baryon mass m^ serves as a "natural" subtraction 
scale in BCPT using the infrared regularization. We follow this convention and set this 
subtraction scale equal to m N . 



III. AXIAL CURRENTS 



In writing down the octet axial currents, we follow standard conventions and notations. 
The most general meson-baryon Lagrangian at lowest order is 
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One may obtain vector and axial vector current operators from £ by including vector 
and axial vector sources in the covariant derivatives. The leading O(p ) operator contains 
only baryon fields and the LEC's D and F . Axial currents involving both baryons and 
mesons first appear at 0(p). Additional purely baryonic axial currents appear at 0(p 2 ). 
They arise from the SU(3) SB Lagrangian 
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(11) 



where 



X+ = ^ + x£ + + £x + (12) 



X = (13) 




The LECs d\-^ are expected to be order of unity in our normalization. There are two other 
terms involving \: 

Tr (B^ lh [A^B\) Tr( X+ ) , Tr (Bj^ 5 {A^ B}) Tr( X+ ). (14) 

These terms do not break SU(3) symmetry and can be absorbed into the definition of D, F 
terms in Eq. (||). 

Using £ ,i one obtains the axial current: 



J A = \D Tr (Bw 5 {er A tf + ?T A i, B}] 

+ \F Tr (B^^T^ + ^T A i, B] 

+ \d^ Tr (B-f7s{ZT A ? + ^^ X+}b) 

+ \d 2 ^ Tr (Bj^T^ + gT A £)B X +) 

+ \d^ Tr (Bn 5 X + BUT A e + ?T A 0) 

+ \d^- Tr (Br l5 B{£T A e + ^T A £, X+ }) 

+ \ Tr (b^T*? -gT A £,B]) 

+ fi^Tr T A ((<9 M E)t£ - <9 M ££+) . (15) 

Renormalized matrix elements of J A between octet baryon states may be written as 
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where the first term on the right hand side is the lowest order one. In terms of the D and 
F coefficients, this term reads for different octet states 

uT = (D + F), 

1-N2 _ _2_ D 

o&i- = (D-F), 



1 
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V 
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oi^ = ~(D + -3F) 



a\V* = --^(D-3F) 



<V- =(D-F), 
«Ss- = ^ D + F ^ = ^4+1° (18) 

where the superscript denotes the corresponding SU(3) indices of the current. 

The second term arises from the SB terms in Eq. ([11]). The coefficients a^ are 
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where 



*S£ = 4(d + /) = >/2aggo (19) 



d=*±*. / = ^i. (20) 



The remaining terms arise from the loops of Figs 1 and 2. The coefficients A^, f3^ , 
7y and 6*^ are given in Tables 0-[ITl] and Eq. fl2"8|) below. In presenting the loop results, 
we give complete expressions for the infrared integrals /. When fitting the LEC's D, F, 
and di, . . . ,c?4, however, we include only the pieces occuring through 0(p 3 ). In doing so, 
we follow approach used in Ref. in making the chiral expansion of the loop integrals. The 
denominator is always kept intact while we expand the numerator up to order 0(p 3 ) only 
after finishing the loop integral explicitly. Meanwhile we never expand terms like 

(4"<K) ±1/2 (21) 

in order to preserve the analyticity properties of the integrals. 

Wavef unction renormalization 

The third term in Eq. fll6|) arises from the wave function renormalization in Fig. 1. We 
have 

Zi = 1 + X u I d + X u I d + Ajj I d , 



Z t Z 3 = 1 + ASiJ + AJ/f + \%n (22) 



The coefficients \ x - (X = n, K, rf) are collected in Table |. The function l£ is defined as 

if = \{Amlm\jf (0) - A x - 2m x [I x (m%) + m x J x (0)]} (23) 

where the expressions of J x (0),Ax,I x (m%) are: 

Ax = (^) 2 log/, 2 , (24) 
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where \i = mx/m^. 

Vertex corrections 

The fourth term comes from the vertex correction diagram Fig. 2a. The coeffients /M 
are collected in Table [TT]. The function I x is defined as 

I* = -A. Y - m\lX(ml) + mi Jf - 2^ + -J&y 

where the last two terms arise from expanding the factors -§5f , J]~2 aroun d D — 4 in the 
scalar integrals in the appendix of Ref . @ . 

The fifth term is the vertex correction from the tadpole diagram in Fig. 2b. The 



coefficients 7^ are presented in Table |T|. The function I x = A 



±x- 



Seagull graphs 

The last term in Eq. ( |I6|) arises from the 0{p) one- meson operators in J^. The relevant 
Feynman diagrams are Fig. 2c and 2d. The contribution from these diagrams is entirely of 
recoil order, vanishing in the m^ — > 00 limit. It was not included in the previous HBCPT 
analyses which worked to leading order in the l/m^ expansion. For this contribution, we 
have 

If = -\m x I*{ml) (27) 

*5 = -^ (28) 

with X = K, 77. 

It is straightforward to verify that we recover previous results in Refs. 0fJ] if we use the 
relation rrii = \vp? k and keep only m? K lnm 2 K terms in Eq. flip]). 



IV. NUMERICAL ANALYSIS AND DISCUSSIONS 

From the expressions of I* b c d we find that the full one loop result Eq. fll6|) contains 
terms of 0(p 3 ) through 0(p 5 ). The terms of odd chiral order (p 3 , p 5 ) are entirely non- 
analytic, whereas the loops yield both analytic and non-analytic contributions of even chiral 
order (p 2 , p 4 ). Recoil order contributions first occur at 0(p 3 ). The contributions at this 
order have the form given in (||D with k = 1 and 

f{fjt) = arccos(-|) x (4 - /x 2 ) ±1/2 . (29) 

Although the /(/x) is non-analytic in the complex plane, one may nevertheless expand it in 
powers of /x 2 about /x 2 = along the real axis. When multiplied by the prefactor (mx/A x ) 2 x 
/x of (@), the leading term in the series scales as m x /A 2 mM, making it of chiral order p 3 . The 
remaining terms - corresponding to successively higher orders in p - represent sub-leading, 
one-loop recoil corrections. We note, however, that this series cannot be reproduced by any 
combination of operators in the effective Lagrangian. Each term in the series scales as an 
odd power of m x , that is, as a fractional power of quark mass. Thus, the infinite series of 
recoil corrections given by the factors in Eq. fl29|) is unambigously associated with loops. 
In HBCPT, one would compute these sub-leading corrections order by order in 1/mjv and 
would be forced to truncate at some order. 

Loop contributions involving even powers of p cannot be disentangled from terms in the 
effective Lagrangian. For example, both the wavefunction renormalization diagrams and the 
vertex corrections generate analytic contributions of 0(p 2 ). Since these contributions are 
quadratic in m x , one could just as well absorb them into the SB terms of Eq. ( pi]) . Similarly, 
at 0(p 4 ), one encounters a new set of SB contributions generated by the Lagrangian 

4 

£2 = -rf {d 5 Tr (5>f>y 5 x + A ti x + B) + d 6 Tr (b^A^BxI) 

x 

+d 7 Tr (B^ l5X+ B{A^ X+ }) } + d 8 Tr (b^A^BA^) +■■■ (30) 

where we have included only a few of the relevant terms. Since the number of LEC's 
appearing at this order is larger than the number of available data, we truncate at 0(p 3 ) in 
order to avoid introducing model assumptions. 

To this order, we determine the LEC's D, F and di, . . . , dj from hyperon semileptonic 
decay data ||, presented in terms of axial vector couplings in Table |V[ Since m u ~ m d « 
m s , we set m n = in performing our numerical fits. We also follow Ref. [Q and enhance the 
errors by 0.2 to avoid the biasing the fit to the precisely known D + F value from neutron 
beta decay. The tree level best fit yields D = 0.78, F = 0.47 and F/D = 0.60 with a 
X 2 = 0.1 for six data points as presented in Table [TV]. In obtaining fits at 0(p 2 ) and beyond, 
we follow two different procedures for treating the analytic loop terms: scheme B, in which 
all the 0{p 2 ) analytic loop contributions are kept explicitly, as in Refs. [0,0 ; and scheme C 
in which these analytic terms are effectively absorbed into the 0(p 2 ) LEC's. 

It is interesting first to truncate at 0(p 2 ) and determine the impact of the SB contri- 
butions. In Ref. 0, where these terms were omitted, the best fit values for the LEC's are 
D = 0.56 and F = 0.33 with F/D = 0.6. Inclusion of the SB terms shifts these values to 
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D = 0.55 and F = 0.41 (D = 0.51, F = 0.37) in scheme B (C), a 10-25% shift. Similarly, 
in when the recoil corrections are included but SB terms omitted, the values for D and F 
are both reduced by roughly 30% from results in Ref. (see Tables [V] and |V1J| ). The full 
0(p 3 ) results yield values of D and F nearly 25% smaller, with F/D remaining close to 0.6. 
The dominant effect arises from inclusion of recoil. 

It is also instructive to determine the numerical importance of including the full recoil 
factors appearing in Eq. fl2"9p. To that end, we make the replacements 

arccosf — — ) — >■ — (31) 






2/ 2 



- n 2 -> 2 (32) 

which amounts to retaining only the leading 1/mjv corrections appearing at 0(p 3 ). Taking 
this limit is equivalent to working to first order in 1/tojv with HBCPT. In this case, the best 
fit values for D and F are essentially unchanged from the 0(p 3 ) fit for scheme C, while the 
SB LEC's dx, . . . , d& shift somewhat. Since the impact of the SB terms relative to the recoil 
corrections is small, it appears that retention of the leading recoil corrections is sufficient to 
determine the convergence behavior of the expansion. 

More significantly, the 0(p 3 ) contributions, arising entirely from recoil effects, are as 
large if not larger than the 0(p 2 ) contributions. The relative importance of each order is 
shown in Tables [VI| and |VH| . To illustrate, we write here the results for a few representative 
channels (in scheme C): 

#4 = 0.658[1 + 0.419 + 0.495] = 1.26 (33) 

g £ A = -0.488[1 - 0.252 + 1.07] = -0.88 , (34) 

where the terms inside the square brackets represent the relative size of the order p°, p 2 , and 
p 3 contributions, respectively. A similar pattern holds for the other octet axial vector matrix 
elements. Far from improving the convergence behavior of the octet-only chiral expansion 
of the axial currents, inclusion of recoil corrections makes it worse. 

In order to explore further why the chiral corrections through 0(p 2 ) are so significant, 
we collect the values of loop integral functions in Eq. (|i6|) in Table [IX|. First, we note 
that the contribution of the vertex correction from Fig. 2a is suppressed due to its small 
coefficients (3$, which are cubic functions of D, F. Although the coefficients XijC^j are also 
cubic in D,F, the coefficients of A^ are big as can clearly seen in Table |. Consequently, 
wavefunction renormalization has a significant impact. Moreover, the contributions from the 
self-energy, tadpole, and seagull diagrams all have the same sign as the the tree level axial 
couplings. These contributions add constructively. In addition, the coefficients of tadpole 
diagram and seagull diagrams are linear function of D, F, so they are enhanced relative to 
the other loops in this respect. 

The relative size of the recoil corrections requires further explanation. To illustrate, 
consider the seagull contributions. Naively, the latter ought to be suppressed by roughly 
"W^tv ~ 1/2 relative to the 0{p 2 ) loop effects. However, the presence of the arccos(— /x/2) 
in these loops generates an additional numerical factor of tt for these contributions at leading 
order in the 1/mjy expansion. It is both the large size of the kaon mass and this numerical 
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factor which are responsible for the large size of the 0(p 3 ) effects. Moreover, such numerical 
enhancement factors appear at higher orders as well. For example, the 0(p 5 ) contributions 
generated by wavefunction renormalization are also proportional to arccos(— /x/2). Thus, we 
would expect the pattern shown in Eq. (|33|) to persist to higher ordersp 1 ]. 

Finally, we illustrate the practical consequences of axial vector non-convergence by con- 
sidering the strange quark contribution to the nucleon's spin, As. As shown in Ref. ||11|| , 
one may express As in terms of the polarized structure function integrals 

r p , n = dx gi' n (x) (35) 

Jo 

as 

As= 3 -[T p + T n }-^gi (36) 

where g$ is the axial vector coupling associated with the matrix element {p\J^\p). The 
combinations of LEC's required for this matrix element are 

a 8 pp = ^=(3F-D) (37) 

/%* = -^(|D» - 2D*F) (38) 

^ = 2l73 (3F " D)3 (39) 

< = ^=(^ 2 -2d 4 ) (40) 



t P K = ~. 7^ = (41) 



0% K ' V = -*&*'" ( 42 ) 

Using our results in scheme C, we obtain 

^ = 0.11[1 + 0.55 + 1] (43) 

where the terms in the brackets correspond to the order p°, p 2 , and p 3 contributions. Using 
this result and the world average data for the r pn , we obtain 



* 1 The problem of the large kaon mass in SU(3) CPT with baryons has also been addressed in 
Ref. 0. 
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As = 0.14 -0.16(1 + 0.55 + 1] (44) 

where we have omitted the experimental error bars in the first term taken from the polarized 
deep inelastic scattering (DIS) data. The second term represents the contribution from 
g%, broken by successive orders as above. We do not quote a total for As given that 
the 0(p 3 ) contribution from g% is as large as both the O(p ) term as well as the first 
term on the RHS of Eq. (@). Given the poor convergence behavior of the expansion 
of $f , extraction of As from polarized DIS data is problematic. In contrast, extractions 
of As from semi-inclusive measurements performed by the Hermes collaboration or elastic 
neutrino-nucleon scattering are not plagued by large SU(3)-breaking uncertainties. Whether 
inclusion of decuplet intermediate states reduces these SU(3)-breaking uncertainties requires 
further study. 
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Figure Captions 

FIG 1. Feyman diagrams for the wave function renormalization. The dashed and solid line 
denotes the pseudoscalar meson and baryon respectively. 

FIG 2. The loop diagrams for the chiral corrections to the axial charge. The filled circle is 

the insertion of the axial current in Eq. flTEp. 
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TABLES 





kaon loop 


r/ loop 


7T loop 


Apn 


fD 2 - ADF + 6E 2 


±L> 2 - 2DF + 3F 2 


3(L> + F) 2 


Aae- 


§L> 2 + 8E 2 


4^2 

3-^ 


| D 2 + 4F 2 


A s o H - 


fD 2 + ADF + 6F 2 


ID 2 + 2DE + 3F 2 


3(D - F) 2 


ApA 


ID 2 - 2DF + 9F 2 


Id 2 -df+\f 2 


ID 2 + 3DF + §E 2 


A AS - 


ID 2 + 2DE + 9F 2 


Id 2 + df + \f 2 


|£> 2 - 3DE + §E 2 


A n S- 


^D 2 -2DF + 5F 2 


Id 2 -DF + §E 2 


^■D 2 + 3DF + ^F 2 


As°h- 


^D 2 + 2DF + 5E 2 


\D 2 + DF+ §E 2 


fD 2 -3DF + ^F 2 


App 


^L> 2 - ADF + 6E 2 


\D 2 - 2DF + 3E 2 


3{D + F) 2 


Aaa 


§L> 2 + 12F 2 


±D 2 
3^ 


AD 2 


Aes 


AD 2 + 4F 2 


3^ 


\D 2 + 8F 2 


Ahs 


f D 2 + 4DE + 6E 2 


±L> 2 + 2DE + 3F 2 


3{D - F) 2 



TABLE I. 



The coefficients X^ for the wave function renormalization. 





kaon loop 


7] loop 


7T loop 


Ppn 


-D*+D Z F-3DF Z +3F A 
3 


D d -5D^F+3DF^+9F- i 
12 


\(D + Ff 


/3as- 


-Te^ + Te^ 2 


\/6 n 3 
9 u 


^(f - 2F 2 ) 


fep~- 


D' 1 +1)^ F+3Di^ +3F d 
3 


D i +hD' z F+3DF' L -9F i 

12 


\{D-Ff 


PpA 


-5D 6 +15D'^F+9DF'^-27F- i ] 


^[-ID3 + |^ 2 ] 


&D{D 2 - F 2 ) 


0AZ- 


-5D d -15D^F+9DF^+27F d ] 
6VE 


^[-^ + IDF 2 ] 


^D{D 2 - F 2 ) 


PnT,- 


D- i +D^F+3DF i +3F i 
6 


-\D Z + \D 2 F - \DF 2 


D J -2D' z F+3DF' 2 +6F- i 
6 


/W 


-D- i +D- i F-3DF z +3F J 

6\/2 


D A +AD^F+3DF-> 
QV2 


D i +2D z F+3DF z -&F- i 
%\[2 



TABLE II. The coefficients /3*j for the vertex correction of Fig. 2a. 
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kaon loop 


77 loop 


■K loop 


7pn 


l 

2 





-1 


7ae- 


1 
2 





-1 


7=0 =- 


1 
2 





-1 


7pA 


3 

4 


3 

8 


3 

8 


7ah- 


3 

4 


3 

8 


3 

8 


7„s- 


3 

4 


3 

8 


3 

8 


7s»H- 


3 

4 


3 

8 


3 

8 



TABLE III. 



The coefficients 7/7 



for the tadpole diagram Fig. 2b. 





Experimental data 


tree level fit 


q A 

Hpn 


(1.2573 ±0.0028) 


1.253 


Sae- 


v /§(0.742± 0.018) 


0.64 


aios- 1 


- 


0.31 


a A 


~y§(0.718± 0.015) 


-0.90 


a A 


yK0.25i0.05) 


0.26 


9ny,- 


(0.340 ±0.017) 


0.31 


5s«H- 


^(1.278 ±0.158) 


0.89 


D 


- 


0.78 


F 


- 


0.47 


F/D 


- 


0.60 


x 2 


— 


0.1 



TABLE IV. Experimental data and O(p ) fits for the axial charge from hyperon semileptonix 
decays. The value in the bracket is the experimental value of gi/fi- The channel with ' is the 
prediction. 



15 





Data 


Tree level fit 


One loop 0(p 2 ) Fit B 


One loop 0(p 3 ) Fit B 


q A 

Hpn 


(1.2573 ±0.0028) 


1.253 


1.245 


1.25 


Sae- 


v /§(0.742± 0.018) 


0.64 


0.60 


0.62 


4. S - f 


- 


0.31 


0.126 


0.29 


9pA 


-^/§(0.718± 0.015) 


-0.90 


-0.90 


-0.89 


a A 
9as- 


yK0.25i0.05) 


0.26 


0.31 


0.30 


9^- 


(0.340 ±0.017) 


0.31 


0.35 


0.34 


5s«H- 


^(1.278 ±0.158) 


0.89 


0.905 


0.90 


D 


- 


0.78 


0.55 


0.41 


F 


- 


0.47 


0.41 


0.26 


F/D 


- 


0.60 


0.75 


0.63 


x 2 


- 


0.1 


0.01 


0.022 


di 


- 


- 


-1.08 


-2.75 


d 2 


- 


- 


0.505 


0.88 


ds 


- 


- 


-0.574 


-0.65 


di 


— 


— 


0.82 


-0.11 



TABLE V. Our fit with Scheme B up to 0(p 2 ), 0{p 3 ). The channel with ' is the prediction. 





Full fit results 


Tree level only 


0{p 2 ) only 


0(p 3 ) only 


Q A 

Hpn 


1.25 


0.671 


0.245 


0.334 


9ir,- 


0.62 


0.33 


0.079 


0.211 


riSL-t 


0.29 


0.143 


-0.031 


0.178 


9 P A 


-0.89 


-0.489 


0.123 


-0.524 


9as- 


0.30 


0.157 


-0.063 


0.206 


9nX~ 


0.34 


0.143 


0.053 


0.144 


SS«H- 


0.90 


0.474 


-0.158 


0.584 



TABLE VI. The separation of our full up to 0(p 3 ) fit results with Scheme B into tree level, 
pure 0(p 2 ), and 0(p 3 ) pieces for the sake of the discussion of convergence of the chiral expansion. 
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Data 


Tree level fit 


One loop 0(p 2 ) Fit C 


One loop 0(p 3 ) Fit C 


q A 

Hpn 


(1.2573 ± 0.0028) 


1.253 


1.26 


1.26 


Sae- 


v /|(0.742± 0.018) 


0.64 


0.64 


0.61 


4. S - f 


- 


0.31 


0.24 


0.22 


9pA 


-y§(0.718± 0.015) 


-0.90 


-0.85 


-0.88 


a A 
9as- 


yK0.25i0.05) 


0.26 


0.34 


0.30 


9^- 


(0.340 ± 0.017) 


0.31 


0.35 


0.34 


5s«H- 


^(1.278 ±0.158) 


0.89 


0.88 


0.90 


D 


- 


0.78 


0.513 


0.39 


F 


- 


0.47 


0.370 


0.26 


F/D 


- 


0.60 


0.72 


0.67 


x 2 


- 


0.1 


0.004 


0.013 


di 


- 


- 


-1.17 


-2.73 


d 2 


- 


- 


0.60 


0.82 


ds 


- 


- 


-0.62 


-0.54 


di 


— 


— 


0.84 


0.094 



TABLE VII. Our fit with Scheme C up to 0{p 2 ), 0(p 3 ). The channel with ' is the prediction. 





Full fit results 


Tree level only 


0{p 2 ) only 


0(p 3 ) only 


Q A 

Hpn 


1.26 


0.658 


0.276 


0.326 


9ir,- 


0.61 


0.318 


0.095 


0.197 


riSL-t 


0.22 


0.12 


-0.042 


0.142 


9 P A 


-0.88 


-0.488 


0.136 


-0.528 


9as- 


0.30 


0.17 


-0.089 


0.219 


9nX~ 


0.34 


0.12 


0.101 


0.119 


SS«H- 


0.90 


0.465 


-0.135 


0.57 



TABLE VIII. The separation of our full up to 0(p 3 ) fit results with Scheme C into tree level, 
pure 0(p 2 ), and 0(p 3 ) pieces. 
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kaon loop 


rj loop 


la 


0.21 


0.27 


h 


-0.23 


-0.237 


Ic 


0.167 


0.230 


Id 


0.34 


0.424 


Ja 


0.533 


0.504 


A 


-0.23 


-0.237 


Urn},) 


-1.37 


-1.53 



TABLE IX. The values of loop integral functions in Eq. (|lC|). 



~e 



y 



-o 



